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ABSTRACT 



Context. Metal-poor globular clusters (GCs) can provide a probe of the earliest epoch of star formation in the Universe, being 
I | the oldest stellar systems observable. In addition, young and intermediate-age low-metallicity GCs are present in external 
, galaxies. Nevertheless, inferring their evolutionary status by using integrated properties may suffer from large intrinsic 
' uncertainty caused by the discrete nature of stars in stellar systems, especially in the case of faint objects. 
\Q \ Aims. In this paper, we evaluate the intrinsic uncertainty (due to statistical effects) affecting the integrated colours and 
f^} . mass-to-light ratios as a function of the cluster integrated visual magnitude {My 1 ), which represents a quantity directly 
measured. We investigate the case of metal-poor single-burst stellar populations with age from a few million years to a likely 
upper value for the Galactic globular cluster ages (~15 Gyr). 

Methods. Our approach is based on Monte Carlo techniques for randomly generating stars distributed according to the cluster's 
O ' mass function. 

f~i Results. Integrated colours and mass-to-light ratios in different photometric bands are checked to be in good agreement with 
the observational values of low-metallicity Galactic clusters; the effect of different assumptions on the Horizontal Branch (HB) 
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morphology is shown to be not relevant, at least for the photometric bands explored here. We present integrated colours and 
mass-to-light ratios as a function of age for different assumptions on the cluster total V magnitude. We find that the intrinsic 
uncertainty cannot be neglected. In particular, in models with My 1 = —4 the broad-band colours show an intrinsic uncertainty 
so high as to prevent precise age evaluation of the cluster. The effects of different assumptions on the initial mass function 



■ and on the minimum mass for which carbon burning is ignited on both integrated colours and mass-to-light ratios are also 
analyzed. Finally, the present predictions are compared with recent results available in the literature, showing in some cases 
non-negligible differences. 
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1. Introduction 

A key issue in astronomy is to determine ages and chemi- 
cal compositions of those stars and stellar systems that 
are needed to reconstruct the formation and evolution 
of galaxies. Toward this goal, the analysis of the stellar 
clusters population in external galaxies and in the Milky 
Way is fundamental for tracing the history of the par- 

Send offprint requests to: Gabriella Raimondo, e-mail: 
raimondoOoa-teramo . inaf . it 

* All the tables (i.e. for models with My* = -4, -6, -8) 
are available in electronic form at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http:/ /cdsweb.u-stra sbg.fr/cg iT)in/qcat?J/A-|-A/| 

** All the figures are available as colored figures in the elec- 
tronic edition of the Journal. 



ent galaxy (see, e.g. ICote et al.lll998t IWest et al.ll2004|) . 

Contrarily to the situation in the Galaxy, in which mas- 
sive (« 10 5 -7- 10 6 stars) stellar clusters are old ( ^ 10 Gyr) 
and metal poor ([Fe/H] ^ —0.7), observations indicate 
that in the Magellanic Clouds, in the Local Group galax- 
ies, and even in galaxies beyond the Local Group, massive 
stellar clusters spread a large range of both metallicity and 
age; thus they are considered the main indi cators of stel- 
lar formation events in the galaxies history ijLarsenl feoOO: 
iMatteucci et al.l l2002t lHarrisI 1200.^ . Hereinafter we will 
use the term "globular clusters" (GCs) to signify massive 
clusters of whatever age and chemical composition. 

GCs have the advantage to be bright objects, and 
then to be easily observed beyond the Local Group. 
In addition, they are thought to be simple stellar pop- 
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ulations (SSPs), consisting of a gravitationally bound 
group of stars born at nearly the same moment, and 
with a nearly identical chemical composition. The in- 
tegrated broad-band colours, line indices, and mass-to- 
light ratios, we observe from those systems, are the 
unique observational tools to understand their proper- 
ties. Since the pioneering work by Tinsley (1972), different 
roups have developed population synthesis mo dels, e.g. 
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therein) in order to interpret such observables. 

Besides the systematic uncertainties due to different 
sets of stellar evolutionary tracks and different spectral 
libraries-used to transform the models from luminosity 
and effective temperature to observable quantities- (see 
e.g. [Chariot et alJ Il99ri iMarastonl Il998t IBrocato et alJ 
120001 iBruzual fc Chariot! l2003t lYil l2003h . broad-band 
colours may suffer from large intrinsic fluctuations caused 
by the discrete nature of the number of stars in the sys- 
tem. The first studies we re carried out in the optical by 
Barbaro fc Bertellil 19771) - for population I clusters - and 
Chiosi et all l|l988(l - fo r intermediate (Z=0.00 1) and so- 



lar metallicity -, while ISantos fe Frogell (|l997t) analyzed 
the case for near infrared (NIR) bands. Among other re- 
sults, the quoted authors concluded that it is necessary to 
include stochastic effects when deriving ages and metal- 
licities from integrated broad-band colours. 

Most studies normalized their theoretical predictions 
to the tot al number of s tars or tota l mass in the clus- 
ter, while IBrocato et all l|l999l l2000|) derived the mean 
broad-band colours and the corresponding dispersions as 
a function of the cluster visual magnitude (My*) for se- 
lected values of ages. This approach directly links theory 
and observations and it is crucial when the cluster age 
and metallicity are inferred from the observed broad-band 
colours, especially for clusters at the faint end of the GC 
luminosity function. 

In this paper we extend the investi gations by 
Bro cato et al 1 lll999l) and lBrocato et akl l)2000|) and we an- 
alyze stochastic effects not only on broad-band colours, 
but also on mass-to-light ratios as a function of the 
adopted cluster visual magnitude, My*. We provide a 
comprehensive study on this problem by investigating the 
time-evolution of both integrated colours and mass-to- 
light ratios for a fine grid of stellar ages and for three 
different values for My 1 . The analysis is carried out using 
new single-burst low-metallicity models (,£ =0.0002) based 
on the updated stellar models database bv ICariulo et al.1 
( 200^). We choose this metallicity because the metal-poor 
GCs can provide a probe of the earliest epoch of star for- 
mation in the Universe, being the oldest stellar systems 
observable. In addition, young and intermediate age low 
metallicity GCs are present in external galaxies (Larsen & 
Richtler 1999, 2000). For these reasons, and also because 
an analysis of young metal-poor clusters gives an idea on 



how old GCs appeared when they formed, we explore a 
wide range of ages (7.7 < log [age (yr)] < 10.3). We note 
that we did not take into account redshift effects thus our 
calculations can be used only for objects with a redshift 
lower than about 0.1. 

The results are compared with a sample of low- 
metallicity clusters in the Galaxy with different HB mor- 
phology, chosen as prototypes of the old stellar popula- 
tions studied in this work. The present theoretical predic- 
tions are also compared with recent results available in the 
literature showing in some cases non-negligible differences. 

We also discuss the influence of the adopted Initial 
Mass Function (IMF) on integrated colours and mass-to- 
light ratios and the effect of changing the maximum mass 
(M up ) for which carbon burning is not ignited, due to ef- 
fects of the degenerate pressure and neutrino energy losses 
in the core. The assumption of a fixed My* can lead to 
a peculiar behaviour when varying the shape of the IMF, 
since an adjustment of the total number of stars might be 
required to keep the My' value fixed. 

The layout of the paper is the following. In Section [21 
the ingredients of the stellar population synthesis code are 
outlined, together with a brief description of the method 
adopted to derive the integrated quantities. In Section |31 
we show the comparison of the theoretical results with 
selected observations of galactic globular clusters. Then, 
we discuss the uncertainties affecting integrated colours 
(Section 0} and mass-to-light ratios (Section |SJ), together 
with a comparison with previous works. 



2. Description of the code 

Synthetic CMDs and magnitudes presented in this paper 
are base d on the stellar population synt hesis code devel- 
oped hv IBrocato et~aTl ill 9991 l200(t . and iRa.imondo et alJ 



(2005) 1 . In this section, we briefly describe the main in- 
gredients and recall the method used to derive integrated 
magnitudes and colours, referring to the cited papers for 
more details. 



2.1. Ingredients 

The present SSP models rely on the evolutionary tracks of 
the "Pisa Evolution ary library" 2 for masses M > 0.6M Q 
l(Cariulo et aLlEooih . The input p hysics adopted in th e 
models has already been discussed in lCariulo et alJ < 2004h . 
We only point out here that the models take into ac- 
count atomic diffusion, including the effects of gravita- 
tional settling, and thermal diffusi on with diffusion coeffi- 
cients gi ven by iThoul et al J (1994); radiative ac celeration 
(see e.g. iRicher et al.lll998t iBichard et al. J |2002h is not in- 
cluded. The effects of rot ation (see e.g. iMaeder fc ZahrJ 
Il998t IPalacios et al"] l2003) arc also not included. 



http:/ /www.oa-teramo.inaf.it/SPoT| 



http://astro.df.unipi.it/SAA/PEL/ZO.html Data files are 
also available at the CDS. 
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Convective regions, identified following the 
Schwarzschild criterion, are treated with the mixing 
length formalism. We use throughout the canonical 
assumption of inefficient overshooting, so the He burning 
structures are calculated according to the prescriptions 
of canonical semiconvection induced by the penetra- 
tion of convective e lements in the radiative region 
i Castellani et aljfl985|) . The efficiency an d presence of a 
mild overshooting are still open questions ijBarmina et alJ 
l2002t iBrocato et all Eooij) : however, as discussed in |YJ 
l)2003|) . a modest amount of overshooting (i.e. Hp ~ 0.2, 
see also Brocato et al. 2003) influences only integrated 
colours with ages 1.5 Gyr for a maximum amount 
of «0.1 mag. The models span the evolutionary phases 
from the main sequence up to C ignition or the onset of 
thermal pulses (TP) in the advanced Asymptotic Giant 
Branch (AGB) in the mass range 0.6-^11 M Q . This allows 
us to calculate stellar population models in the age range 
w 50 Myr 20 Gyr. 

Beyond the early-AGB phase, thermally pulsating 
(TP) stars have been simul ated in the synthesis code us- 
ing th e analytic formulas of IWagenhuber & Groenewegenl 
( 199j|) which describes time evolution of the core mass, 
and luminosity of TP stars. These formulae include three 
important effects: (i) the fact that the first pulses do 
not reach the full amplitude, (ii) the hot bottom burn- 
ing process that occurs in massive stars, and (m) the 
third dredge-up. The effective temperature (T e ) of each 
TP-AGB star has been evaluated using prescriptions by 
iRenzini fc Vol] ilSsl . considering the appropriate slope 
dlog(L/ Lq) / dlog(T e ) of the adopted evolutionary tracks. 
The analytic procedure ends up providing the time evolu- 
tion of the te mperature and lu minosity for a given mass 
(see for details IR.aimondo et ahlEoOSh . 

Mass loss affecting Red Giant Branch (RGB) stars and 
early-AGB star s has been tak en into account following 
prescriptions bv lReimersl fll975|) : 



M R = -4 • lO" 13 ^ • LR/M, 



(1) 



while during the TP phase we adopted the 
iBaud fe Habind l|l983j) mass-loss rate: 



M BH = fiLR/M e . 



(2) 



Here, L, R, M, M c are, respectively, the star luminosity, 
radius, total mass, and envelope mass in solar units; 
H = — 4 ■ 10 _13 (M ej o/M), being M ei0 the envelope mass 
at the first TP. Eq. is a modification of the Reimers 
formula with t]r — 1 which also includes a dependence 
on the actual ma ss envelope. The initial-f inal mass rela- 
tion is applied bv iDominguez et aJl (|l999h . and no white 
dwarfs (WD) more massive tha n 1 . 1 Mq are accepted 
l|Prada Moroni fc Stranieroll2005|) . 

In this paper the TP phase is included but we do not 
adopt any separation between C-rich and O-rich TP-stars. 
All stars are oxygen-rich when they enter the AGB phase. 
Whether or not they become C-stars depends primarily 
on the efficiency of the third dredge-up (TDU) occurring 



on the TP-AGB phas e, and the extent and time-variation 
of th e mass-loss (e.g. IMarigo et al.lll999l IStraniero et alJ 
2003). In low-metallicity stars (Z < 0.004), the amount 
of oxygen in the envelope is so low that a few thermal 
pulses are sufficient to convert an O-rich star into a C-star 
!|Renzini fc Volilll98l|) . In addition, the smaller the metal- 
licit y the smaller the minim um mass for the onset of TDU 
(e.g. IStraniero et alJEool . On the other hand, TP-AGB 
stars may experience episode of strong mass loss, that in 
the case of low mass stars may cause a reduction of the 
envelope mass that may delay or even preven t the TDU oc- 
curre nce and the formation of C-rich stars 1 Marigo et alJ 
1999). In conclusion, the presence of C-rich stars may af- 
fect NIR-bands luminosity and its uncertainty in the case 
of l ow-metallicity, i ntermediate-age massive clusters (see 
e.g. lMarastonlll99'8|l . while at the typical age of Galactic 
globular clusters their presence become more uncertain, 
as confirmed by the fact that AGB stars in GGC are all 
observed to be oxygen-rich, so that carbon does not ap- 
pear to have been dredge-up into the en velops during ther- 
mal pulses l|Lattanzio fc Woo'dl l2003|) . In addition, our 
assumption is also expected to have marginal effect on 
integrated quantities of faint populations, as bright TP- 
AGB stars are statistically less frequent, or even absent. 
Finally, it will be shown that the details and the treat- 
ment of the physical processes at work on the TP-AGB 
phase, as well as their impact on synthetic colours, is still 
uncertain (Section 0J. 

The adopted colour transfor mations f or the standard 
UBVRIJHK band s are from ICastellil (|l999|) . see also 
ICastelli et al ] <1997h . To calculate integrated colours in 
the Hubble Space Telescope (HST) bands (WFPC2 and 
NICMO S systems) we adopte d the colour tr ansformations 
bv lOriglia fc Leithererl l)2000h based on the iBessell et alJ 
l|l998h stellar atmospheric models. 

As well known the mixing length parameter (a) gov- 
erns the efficiency of convection in the convective envelope 
of a stellar structure. The a parameter used for calculat- 
ing the evolutionary tracks adopted in the present work 
has been calibrated in such a way that the isochrones re- 
produce, with the adopted colour transformations, the ob- 
served RG branch colour of GCs with the proper metal- 
licity and age. This is evident from Fig. [5] in which our 
synthetic models nicely reproduce the RGB colours of the 
selected GCs. The a parameter needed in the present evo- 
lutionary models to obtain this agreement is a «2.0, but 
it's worth noticing that a is a free parameter, sensitively 
dependent on the chosen color transformations and on 
all the adopted physical inputs which affects the effective 
temperature of a model. ICariulo et al.l 1 20041) showed that 
tracks with the same physical assumptions for metallicity 
up to Z«0.001 reproduces the RGB colours of galactic 
GCs for the same a values, while metal-rich (Z>0.001) 
and standard solar models require an a s lightly lower (see 
ICiacio et akllT997tlCa,stella.ni et aljEiol. 

The ver y low mass tr a cks (V LM, M < 0.6 Mq) are 
taken from iBaraffe et all l|l997|) . The tracks have been 
already transformed by the authors to the observational 
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plane in the Johns on-Cousins sy s tem a dopting the colour 
transformations bv lAllard et al . (1997). particularly suit- 
able for low mass stars. We checked that low mass mod- 
els satisfactory match, in all the available colours, the 
higher mass mode ls (Fagiolini 20 4). Ho wever, as already 
discussed by e.g. iBrocato et all ((2000^ . VLM stars do 
not contribute to the photometric indices, although their 
contribution to the cluster mass is fundamental. Masses 
lower than the mini mum mass for the c entral H ignition 
~ 0.08 Mq, (see e.g. iBar affe et al]ll997h have been esti- 
mated do not signific antly contribute to the total mass 
of the cluster (see e.g. IChabrier fc Meralll997l) . and thus 
they are not taken into account. 

Post-AGB evolution, until the entrance in the white 
dwarfs (WD) cooling sequence, are not considered be- 
cause the evolutionary t ime is too short to have a sig- 
nificant influence (see e.g lB locker et al.lll997f) . WDs have 
be en included in the code using evolutionary models 
bv lSalaris et"all §000). They have been transformed by 
adopt ing the atmospheric models bv lSaumon fc .Tacobsonl 
(1999) for DA WDs, which include the treatment of the 
collision-induced absorption of H2 molecules, for T < 
4000 if. For hi g her te mperatures the transformations of 
iBergeron et alJ l[t995) have been used. As we will dis- 
cuss in t he following our standa rd model adopts 6.5 Mq 
as M up , l|Dominguez et alJll999h. A s expe cted, and as al- 
ready noted by e.g. lAngeletti et al the contribu- 
tion of the WD population to optical and NIR photometric 
indices is negligible although they significantly contribute 
to the total mass of the cluster. 

The IMF of lKroupal l)2002l) is adopted in the mass in- 
terval 0.1 < M/Mq < 11, unless explicitly stated oth- 
erwise. To simulate the mass distribution of stars in the 
synthetic CMD we use a Monte Carlo method: the position 
of each randomly created star in the logL/L© vs. logT e 
diagram is determined for each given age. As already dis- 
cussed, the chemical composition is fixed to Z = 0.0002 
and Y = 0.23. 

Fig-fflshows, as an example, synthetic CMDs (without 
simulation of photometric errors) for the selected chemical 
composition, a total absolute visual magnitude My* = 
—6 mag, and three different ages (500 Myr, 2 Gyr, and 
11 Gyr). All the evolutionary phases described above are 
clearly visible. 

2.2. Integrated magnitudes and colours 

To compute integrated fluxes and magnitudes we assume 
that the integrated light from the stellar population is 
dominated by light emitted by its stellar component. This 
implies that i) no source of non-thermal emission are at 
work, ii) the thermal emission by interstellar gas gives no 
sizeable contribution to the integrated flux, and iii) the 
absorption from dust and gas are negligible. On this ba- 
sis, the integrated flux in each photometric filter mainly 
depend on two quantities. The first is the flux fi emitted 
by the i-th star of mass M, age t (stellar system age), Lu- 



is 




(v-D 



Fig. 1. Synthetic CMDs for Z = 0.0002, Y = 0.23 and 
three different ages: 500 Myr (black), 2 Gyr (red), and 
11 Gyr (blue). For each population My* = —6 mag is 
adopted. The simulation of photometric errors is not in- 
cluded. 



minosity I, effective temperature T e //, and chemical com- 
position (Y, Z): 



fi[l(M, t, Y, Z), T eff (M, t, Y, Z), Y, Z] 



(3) 



fi is defined by the stellar evolutionary-tracks library, and 
by the adopted temperature-colours transformation ta- 
bles. 

The second quantity is <L>(M, N) which describes the 
number of stars with mass M in a population globally con- 
stituted of N stars with age t and chemical composition 
(Y, Z). It is strictly related to the IMF. 

A fundamental point of our code is that the mass of 
each generated star is obtained by using Monte Carlo 
techniques, while the mass distribution is ruled by the 
IMF. The mass of each star is generated randomly, and its 
proper evolutionary line is computed by interpolating the 
available tracks in the mass grid. This method is crucial for 
poorly populated stellar systems, as we shall study in the 
following. It ensures to treat undersampled evolutionary 
phases properly; for instance NIR colou rs may be dom- 
inated by a handful of red giant stars ij jjjS anfos Frogell 
1997t IBrocato et alJll999t ICervifio fc Valls-Gabaudll2003l 
Raimondo et al.ll2005|) . ~ 
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In each model, stars are added until to reach a given 
value of the absolute visual magnitude, My, at any age. 
The mass values are generated randomly and distributed 
according to the chosen IMF. It is relevant to note that 
the random extraction of masses is fully independent from 
model to model even if the same input quantities (My, t, 
N, Y, Z, IMF, etc.) are assumed. 

Both the previous quantities (fi and $(Af, N)) are 
combined and integrated by the stellar population code 
to derive the total integrated flux F in a given photomet- 
ric band: 



(4) 



F[N,t,Y,Z] = J2fi, 

i—l 

and the corresponding magnitudes. 

3. Comparison with Galactic Globular Clusters 

Before studying the general behaviour of colours and 
mass-to-light ratios as a function of My, the natural test 
for SSP models is to compare their predictions to the ob- 
served properties of Galactic star clusters. Three old galac- 
tic globular clusters (GGCs) have been selected for check- 
ing our capability to reproduce both the CMD morphology 
of each cluster and the integrated magnitudes: NGC4590 
(M68), NGC7078 (M15) and NGC7099 (M30). They are 
all metal-poor, span a relativity wide range of the vi- 
sual magnitude (Tabled, and show different Horizontal 
Branch morphologies. This allows us to check the effect 
of HB morphologies on the integrated colours in the se- 
lected bands. The [Fe/H] values estimated for the three 
clusters are, respectively, [Fe/H] = -2.06,-2.26,-2.12 
(iHarrislllQQ^ 3 . while the a-elements enhancement can be 
evaluated as \a/F e] an 0.3 (see, e.g, the discussion in 
iFerraro et alJfl999|) : thus we can assume Z an 2 ■ 10~ 4 . 
We include M15 due to its brightness (My ~ —9), even if 
it is recognized as a cluster which should have undergone 
a gravothermal catastrophe. This results in a contraction 
of the cluster core while the external regions expand, with 
some stars escape the system. As a possible consequence 
of such a dynamical evolution, observational evidence of 
mass se gregation, and of radial color gradients have been 
found feailyn et all Il989t iDe Marchi fc Paresce I Il994t 
IStetsonlll994h ~ 

Fig. [5] shows the observed CMDs for the selected clus- 
ters. B and V photometry come s from the HST-Snapshot 
Catalog by IPiotto et all ll2002l) . except for M68 whose 
data are from lWalkerl l)l994|) . The synthetic CMD which 
better reproduces the observational data is over-plotted as 
black dots in each panel. For each cluster an age of 11 Gyr 
is assumed; we are not interested in a detailed calibration 
of the cluster ages which is far from our purpose. Results 
for each comparison between synthetic and observed clus- 
ter are, briefly: 
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3 www.physics.mcmaster.ca/resources/globular.html. 
February 2003 version. 



Fig. 2. Comparison between present synthetic (red dots), 
without the inclusion of simulated photometric errors, and 
the observed CMDs of the globular clusters M68, Ml 5 and 
M30 (black dots). The assumed age is 11 Gyr, and the 
adopted chemical composition is Z = 0.0002 and Y — 0.23 
(see text). The values of the distance modulus and the 
reddening obtained from the analysis are also indicated. 



- M68 (NGC4590): We find (m - M) v = 15.29 
and Eb-v — 0.05 in agreement, within the uncertain- 
ties, with the curren t distance modulus and reddening 
determinations (e.g. iHarrid Il996t ICarretta et alJ l2000t 
iDi Criscienzo et al.ll2004(l . The HB morphology is repro- 
duced using rjji = 0.30 with a dispersion a v = 0.08. Due 
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Table 1. Observed and theoretical integrated colours for the selected globular clusters. Values from the Harris' 
catalogue are dereddened according to Reed et al. (1988). The J — K near-IR colour is from Brocato et al. (1990), 
dereddened according to Cardelli et al. (1989). M./Ly is from the compilation by Pryor & Meylan (1993). 



Identity 




U-B 


B-V V-R 


V-I 


J-K 


M/L v 


Observations 


M68 
M15 
M30 


-7.35 
-9.17 
-7.43 


-0.01 ± 0.01 
-0.03 ± 0.02 
-0.00 ± 0.02 


0.58 ± 0.03 0.43 
0.58 ± 0.01 

0.57 ± 0.04 0.39 ± 0.02 


0.88 ± 0.01 

0.73 
0.82 ± 0.04 


0.62 
0.54 


1.6 
2.2 
2.5 


Models 


M68 
M15 
M30 


-7.34 ± 0.05 
-9.18 ± 0.02 
-7.43 ± 0.05 


0.00 ± 0.02 
0.00 ± 0.01 
-0.01 ± 0.01 


0.60 ± 0.02 0.42 ± 0.01 
0.60 ± 0.01 0.420 ± 0.003 
0.59 ± 0.02 0.41 ± 0.01 


0.87 ± 0.02 
0.89 ± 0.01 
0.87 ± 0.02 


0.57 ± 0.07 
0.59 ± 0.02 
0.54 ± 0.05 


1.79 ± 0.07 
1.58 ± 0.03 
1.76 ± 0.07 



to the present uncertainties in the factors that influence 
the HB morphology and in the precise treatment of mass 
loss i n the RG branch (see e.g. iLee et af1ll994t iRev et all 
120011) . the tuning of the adopted value of t]r is just a way 
to obtain the observed HB morphology but it should not 
be interpreted in terms of physical parameters of the stel- 
lar cluster. 

- M15 (NGC7078): The best result is obtained by as- 
suming (m — M)v = 15.50 and Eb-v — 0.09, in agree- 
ment with current deter minations on the cluster distance 
modulus and reddening JHarrisHl996i iDi Criscienzo et alJ 
l2004t iMcNamara et alJl2004h 'The HB morphology is re- 
produced using rjji = 0.44 and cr r) = 0.1. The cluster con- 
tains one of the few planetary nebulae ( PN) known i n 
GGCs: the star K648 identified as a PN bv lPeasel l)l928|) . 
for which Alves et al. (2000) measured an apparent mag- 
nitude of V = 14.73. 

- M30 (NGC7099): Assuming E B -v = 0-05 we obtain 
(m — M)y = 14.93 in agreement with ISandauist et alJ 
( 1999). The HB morphology is r eprod uced using the same 
parameters as in lBrocato et all l)2000j) . i.e. rjn = 0.40, and 
a, n = 0.2. 

Table ^ gives the observed integrated colours for the 
selected clusters, taken from the Catalog of P arameters 
for Milky Way Globular Clusters l|Harris!ll996|) 4 and syn- 
thetic integrated colours obtained from the present models 
by reproducing, within the errors, the observed total vi- 
sual magnitude of each cluste r. Optical observat ional data 
are de-reddened according to lReed et alJ 1 19881) with red- 
dening values taken from Harris (1996). To give an idea of 
uncertainties for the colour s we report the 'residual' values 
computed by iB.eedl lll985l) which result from his homog- 
enization procedure based on measurements of clusters 
colours by various authors. In some cases the availability of 
only one measurement prevents th e determination of a n 
uncertainty. NIR colours are from iBrocato et all l|l990|h 



4 www.physics.mcmaster.ca/resources/globular.html. 
February 2003 version. The integrated colours U — B 
and B — V are on the standard Johnson system, and V — R, 
V — I on the Kron-Cousins system. The values are the simply 
average of results from Peterson (1993) and Reed (1985). 



de-reddened according to ICardelli et alJ l)l989j) with red- 
dening values taken from Harris (1996). 

Theoretical errors correspond to 1 a dispersion evalu- 
ated from 10 independent simulations. Tablenjlists: clus- 
ter NGC number (Col. 1), total absolute V magnitude 
(M*, ot , Col. 2), integrated U-B, B-V, V — R, V-I, 
V - J and V - K colours (C ols. 3-7), and M/L v (Col. 8) 
from lPrvor & Mev lan ( 1993J). The last values are strongly 
model-dependent, as stressed by the authors themselves. 
Hereinafter we will indicate dereddened colours as, e.g., 
U-B instead of (U - B) . 

For each cluster, numerical simulations were performed 
by populating the synthetic CMD until the observed My* 
of the cluster is reproduced (Sect. 2.1). Synthetic colours 
agree with data of all the clusters to within the uncer- 
tainties and theoretical statistical fluctuations. The only 
exception is the V — I colour of M15 that is found to 
be redder than observations. The observed V — I colour of 
this cluster is significantl y bluer than that of other clusters 
with similar metallicity l|Harrisll 19961) . On this regard, we 
recall that the V —I measure available is based on o ne pho- 
toelectric measurement only l|Kron fc Mavalll960j) . and it 
is not possible to estimate the uncertainty (Reed 1985). 
To investigate th i s issu e, we use VI stellar photometry by 
iR.osenberg et alJ l)2000l) . which contains more than 90% of 
the cluster light (My ~ —9.1). By adding the flux of indi- 
vidual stars, we derive an integrated colour V — I ~ 0.85, 
in fair agreement with our p rediction. This value is larger 
than the integrated value bv lKron fc Mavalil l|l960l h even 
if it would be affected by incompleteness at the faint end 
of the observed luminosity function. To deeply investigate 
the origin of the latter inconsistency is well beyond the 
purpose of this paper, what we wish to point out here is 
to notice that the observed integrated V — I colour of such 
a cluster, reported in Tabled could be peculiar and then 
not representative of metal-poor clusters. 

Additional indications may come from the analysis 
of M15 NIR colours. J — K pred iction well agrees with 
the value by IBrocato et all l|l990l) (J - K = 0.67, red- 
dened ) , and is also consistent with data bv lBurstein et all 
( 19 841) who found J — K = 0.62 (reddened). Interestingly, 
iBurstein et alJ l)l984|) also provide the observed (V — K) = 
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2.14 (reddened), that becomes (V — K) = 1.86 if cor- 
rected according to our reddening assumptions. By com- 
paring this value with theory (V — K = 2.05), it appears 
that our models predict V — K colour about 0.2 mag red- 
der than observations. Unfortunately, a similar compari- 
son cann ot be done for M30 an d M68, since they are not 
listed by iBurstein et alJ l(l984h . Instead the authors ob- 
served M92 (NGC6341), whose metallicity is estimated to 
be close to M15, being [Fe/H] = -2.29. Since M92 is 
as bright as My — —8.20, its colours data can be eas- 
ily compared with our predictions as reported in Table [5J 
at age ~ 11 — 13 Gyr. Interestingly, the observed optical- 
NIR colour V — K = 2.13 well agrees with our predic- 
tion, together with all the optical colours U — B = 0.01, 
B — V — 0.61, and V — I = 0.86 (from Harris's catalog). 

In conclusion, except for the peculiarity concerning 
M15, our models are able to properly predict the inte- 
grated colours of let's say 'normal' metal-poor clusters, as 
confirmed by the comparison with observational data of 
M30, M68 (Table P), and the additional cluster M92. 

As noted above, M./L values reported in the last col- 
umn of Table ^ ( upper section), d e rived using isotropic 
King models by IPrvor fc Mevlanl l)l993h . are strongly 
model-dependent and thus the reported values are only in- 
dicative. Mo reover dynamical processes such as star evap- 
oration fe.g. lSpitzerlll987r.lVesperini fc Heggielll997(l . may 
affect the total mass and the inferred IMF shape as a 
function of time. Observed integrated colours of the se- 
lected clusters, as well as theoretical results, indicate that 
they are not influenced by the HB morphology, at least in 
our modeled photometric bands. Shorter wavelength ul- 
traviolet colours are expected to be more sensitive to the 
presence of an extended blue HB. 



i 

V-K 



4. Integrated colours 

In this section we examine the effects on integrated colours 
of stochastic fluctuations for different values of the as- 
sumed total absolute visual magnitude and the effects 
of the still pr esent uncertainty on the IMF shape. The 
iKroupal l|2002f) IMF is adopted for the standard model. 

Results are presented at fixed absolute visual mag- 
nitude to aid comparisons with observational data. For 
these calculations the evaluation of the statistical fluctu- 
ations is fundamental. Statistical fluctuations of broad- 
band colours as well as mass-to-light ratios have been eval- 
uated by computing a series of 10 independent simulations 
for a fixed set of population's parameters (Z, age, IMF, 
...) at fixed My 1 , and assuming a 1 u error. When fluc- 
tuations become large, and the value above is not fully 
representative of colours variations, we extend the anal- 
ysis up to 100 runs. This is especially the case of very 
poorly populated clusters (My = —4) and young ages. 

The standard models are computed assuming: n^ = 
0.3, a v = 0.08 and Kroupa's IMF. As an example, Table 
121 and |3 reports integrated colours as function of age for 
models with My* = —8; Johnson-Cousins colours in Tab|21 
and HST colours in TabEJ TableHlists from left to right: 
age, M^, integrated (U - B), (B — V), {V - R), (V - 
I),{V — J) and (V — K) colours. Table|3]reports from from 
left to right: age, V-F439W2, V-F555W2, V-F8UW2 
for the WFPC2 and V - F110W1, V - F160W1 for the 
NICMOS camera. 



4.1 



Mtf 1 sensitivity 



Fig. 4. V — K colour distribution resulting from 100 sim- 
ulations for a model with 100 Myr, Z = 0.0002, and 
My 31 = -4 (solid) and -8 mag (dotted). 



Our standard synthetic models are calculated for three 
different values of the total magnitude: My* = —8, —6 
and —4. (Note, that all the results for both broad-band 
colours and mass-to-light ratios for = —6, and —4 

models are available as electronic tables at CDS). 

In Fig. |21 the time evolution of selected integrated 
colours are shown for the two extreme cases (My 1 = — 8, 
black circles, and —4, red triangles). In general, broad- 
band colours become redder with age, because of the in- 
creasing fraction of cool stars in the populations. For the 
same reason, the (U — B) colour is quite age-insensitive 
for age 1 Gyr. 

Within the statistical uncertainties, the mean colours 
at different My* do not show significant differences, ex- 
cept in a few cases at young ages, and for optical-NIR 
colours. This is due to the strong variation of the number 
of very bright red stars affecting significantly the total 
luminosity of the cluster. In a SSP with a total bright- 
ness as faint as My 1 = —4, the most of stars occupy 
the MS phase, and only a few (or none) red giants are 
present, including stars burning He in the core and stars 
in the double shells p hase (see for example Fig. 14 in 
iRaimondo et alJ [20051. Hence, the mean V — K colour 
(and to a less extent V — I) is generally bluer than in 
models with My* = —8, where post-MS phase is always 
widely populated. This effect is large in young clusters, for 
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Table 2. Integrated colours in the standard UBVRIJHK photometric filters for our standard model with My * = —8 
and Z=0.0002. 



Age (Myr) 






M^ ot 






U 




B 






B 




V 






V 




R 




V — I 




\ 


- J 




V -K 


50. 




-8. 


D ± 0.1 


-0 


56 


± 





03 


-0 


10 


± 





08 


-0 


02 


± 





08 


0.0 ± 0.2 


0.0 ± 0.3 


0.1 ± 0.5 


100. 




-8. 


] ± 0.1 


-0 


39 


± 





03 





00 


± 





07 





06 


± 





07 


0.2 ± 0.2 


0.4 ± 0.3 


0.7 ± 0.5 


150. 




-8.0 ± 0.1 


-0 


31 


± 





02 





01 


± 





04 





06 


± 





04 


0.2 ± 0.1 


0.4 ± 0.2 


0.7 ± 0.3 


200. 


-8 


00 


± 


0.05 


-0 


22 


± 





01 





09 


± 





05 





14 


± 





05 


0.4 ± 0.1 


0.7 ± 0.2 


1.2 ± 0.3 


300. 


-8 


01 


± 


0.09 


-0 


12 


± 





02 





14 


± 





07 





16 


± 





07 


0.4 ± 0.1 


0.8 ± 0.2 


1.3 ± 0.4 


400. 


-8 


06 


± 


0.08 


-0 


05 


± 





01 





16 


± 





04 





16 


± 





04 


0.41 


± 


0.09 


0.8 ± 0.2 


1.3 ± 0.2 


500. 


-8 


02 


± 


0.05 


-0 


01 


± 





01 





19 


± 





03 





18 


± 





03 


0.44 


± 


0.07 


0.8 ± 0.1 


1.3 ± 0.2 


600. 


-8 


06 


± 


0.08 





00 


± 





01 





26 


± 





03 





22 


± 





03 


0.52 


± 


0.06 


1.0 ± 0.1 


1.4 ± 0.2 


800. 


-8 


03 


± 


0.06 





01 


± 





01 





34 


± 





04 





27 


± 





03 


0.62 


± 


0.07 


1.1 ± 0.1 


1.7 ± 0.2 


1000. 


-7 


98 


± 


0.04 





00 


± 





01 





36 


± 





02 





28 


± 





01 


0.63 


± 


0.03 


1.13 


± 


0.06 


1.6 ± 0.1 


1100. 


-8 


02 


± 


0.06 





00 


± 





01 





39 


± 





03 





30 


± 





02 


0.66 


± 


0.04 


1.16 


± 


0.07 


1.68 


± 0.09 


1500. 


-8 


03 


± 


0.04 





01 


± 





01 





41 


± 





02 





30 


± 





01 


0.67 


± 


0.02 


1.17 


± 


0.04 


1.70 


± 0.07 


1700. 


-8 


00 


± 


0.04 





02 


± 





01 





42 


± 





02 





31 


± 





01 


0.68 


± 


0.03 


1.20 


± 


0.06 


1.73 


± 0.08 


2000. 


-7 


97 


± 


0.03 





03 


± 





01 





44 


± 





01 





32 


± 





01 


0.71 


± 


0.02 


1.24 


± 


0.04 


1.78 


± 0.05 


3000. 


-7 


98 


± 


0.05 





03 


± 





01 





49 


± 





02 





35 


± 





01 


0.75 


± 


0.02 


1.29 


± 


0.03 


1.83 


± 0.05 


4nnn 


o 


04 


± 


0.06 


n 
u 









ni 


n 
u 


JO 




n 
u 


D9 


n 
u 






n 
u 


u ± 


0.79 


± 


0.02 


1.34 


± 


0.03 


1.89 


± 0.05 


5000. 


-8 


01 


± 


0.05 





00 


± 





01 





55 


± 





02 





38 


± 





01 


0.80 


± 


0.02 


1.35 


± 


0.04 


1.89 


± 0.05 


6000. 


-8 


02 


± 


0.03 


-0 


02 


± 





01 





57 


± 





01 





39 


± 





01 


0.83 


± 


0.01 


1.39 


± 


0.02 


1.94 


± 0.03 


7000. 


-8 


03 


± 


0.04 


-0 


02 


± 





01 





60 


± 





01 





41 


± 





01 


0.85 


± 


0.01 


1.43 


± 


0.02 


1.99 


± 0.03 


8000. 


-8 


00 


± 


0.03 


-0 


02 


± 





01 





61 


± 





01 





41 


± 





01 


0.86 


± 


0.01 


1.44 


± 


0.02 


2.00 


± 0.03 


9000. 


-7 


99 


± 


0.02 


-0 


01 


± 





01 





61 


± 





01 





41 


± 





01 


0.87 


± 


0.01 


1.45 


± 


0.02 


2.02 


± 0.02 


10000. 


-7 


94 


± 


0.03 


-0 


00 


± 





01 





61 


± 





01 





42 


± 





01 


0.87 


± 


0.01 


1.47 


± 


0.02 


2.04 


± 0.03 


11000. 


-8 


01 


± 


0.04 





00 


± 





01 





60 


± 





01 





42 


± 





01 


0.87 


± 


0.01 


1.46 


± 


0.02 


2.04 


± 0.02 


12000. 


-7 


99 


± 


0.04 





00 


± 





01 





61 


± 





01 





42 


± 





01 


0.88 


± 


0.01 


1.48 


± 


0.02 


2.06 


± 0.03 


13000. 


-7 


95 


± 


0.05 





01 


± 





01 





62 


± 





02 





43 


± 





01 


0.90 


± 


0.02 


1.51 


± 


0.03 


2.10 


± 0.04 


14000. 


-8 


00 


± 


0.04 





00 


± 





01 





63 


± 





01 





43 


± 





01 


0.90 


± 


0.01 


1.52 


± 


0.02 


2.11 


± 0.03 


15000. 


-8 


01 


± 


0.04 


-0 


02 


± 





01 





62 


± 





02 





43 


± 





01 


0.90 


± 


0.01 


1.51 


± 


0.03 


2.10 


± 0.04 



which colours are more sensitive to t he clu ster luminosity 
[see also Fig. 5 inlSantos fc Frogell l)l997h . and Fig. 8 in 
IBrocato et all l|l999j) ]. 

To clarify this point, Fig. 0] illustrates the V — K 
colour distribution for a SSP model with age 100 Myr; 
for each choice on My 1 , 100 simulations are computed. 
The V — K colour distribution of faint clusters (solid line) 
shows a blue populous peak, due to the large fraction of 
simulations containing mainly blue stars, and a few sparse 
simulations at redder colours containing a few red giants. 
Consequently, the average colour is bluer than that found 
for brighter clusters (Fig. |3J), which show a broadened 
colour distribution shifted towards the red side, since the 
probability to have a large number of red stars is higher 
(dotted line in Fig.^J. We emphasize that in the present 
paper for a given cluster age, the integrated colours at dif- 
ferent cluster brightness are obtained by properly adding 
stars in all evolutionary phases according to the IMF 
and evolutionary time-scales (i.e. keeping fixed the pro- 
portion between the number of stars in different stages) . 
So that, by decreasing their absolute magnitude (mass), 
the natural trend is that the cluster progressively misses 
post-MS stars, and due to the discreteness of stars " if an 
SSP is far less luminous, then there would be no post-MS 
stars and the integrated fluxes would be dominated by up- 
per MS stars with little spre ad in color, which res ults in 
smaller color fluctuations" ( Santos k, Frogell Il997f> . The 



colour behaviou r descr ibed above is in agreeme nt with 
ISantos fc Frogell (|l997h an d IBrocato et alJ lll999h. w hile 
opposite to what found bv IBruzual fc Chariot! l|2003|) . A 

discussion on the nature of such a discrepancy is beyond 
the purpose of this paper, since it would required a deep 
comparison between all the assumptions adopted in all the 
codes (evolutionary tracks, atmospheres, etc.). 

In general, for a further increase of the cluster bright- 
ness, the stochastic effects decrease, the colour distribu- 
tion becomes more regular, and the standard deviation 
squeezes (small spread in colour). 

At a given age, the size of errorbars increases as the 
cluster luminosity decreases; this is due to the small num- 
ber of stars present in faint clusters (My 1 = —4, red sym- 
bols) with respect to bright ones (My* = —8, black sym- 
bols). From the figure it is evident that for all ages the 
statistical errors at My 1 = —4 are so high as to prevent 
precise quantitative evaluations. 

The erro rbars also increase f rom B — V to V — K , 
as shown bv IBrocato et all 1 1999f) . This last finding is re- 
lated again to the expected small number of post-MS stars 
shining in the infrared wavelengths, especially in case of 
low-luminosity clusters (My* = —4). Moreover, in faint 
clusters it may happen that giant stars are not present in 
each of the simulations we computed at fixed age, being 
their appearance highly driven by stochastic phenomena. 
At fixed total magnitude the intrinsic uncertainty lowers 



Fagiolini et al.: Monte Carlo Simulations of Metal Poor Star Clusters 9 
Table 3. As in Table [21 but for WFPC2 (first three columns), and NICMOS1 (last two columns) photometric bands. 



Age (Myr) V - 439W(2) V - F555W(2) V - F814W(2) V - FlWWjl) V - F160W(1) 



50. 





.09 


± 





.08 


0. 


.009 


± 





.005 




0.0 ± 0.2 


0.0 ± 0.3 


0.0 ± 0.5 


100. 


-0 


.01 


± 


0. 


.08 


0. 


.002 


± 





.005 




0.2 ± 0.2 


0.3 ± 0.3 


0.6 ± 0.5 


150. 


-0 


.02 


± 


0. 


.03 


0. 


.000 


± 





.002 




0.2 ± 0.1 


0.3 ± 0.1 


0.6 ± 0.3 


200. 


-0 


.11 


± 


0. 


.05 


-0 


.006 


± 


0. 


.003 




0.3 ± 0.1 


0.6 ± 0.2 


1.1 ± 0.3 


300. 


-0 


.15 


± 


0. 


.07 


-0 


.009 


± 


0. 


.004 




0.4 ± 0.1 


0.6 ± 0.2 


1.2 ± 0.4 


400. 


-0 


.18 


± 


0. 


.05 


-0 


.011 


± 





.003 


0.4 


± 





.09 


0.6 ± 0.1 


1.2 ± 0.2 


500. 


-0 


.21 


± 


0. 


.04 


-0 


.013 


± 





.002 





.42 


± 


0. 


.07 


0.7 ± 0.1 


1.2 ± 0.2 


600. 


-0 


.28 


± 





.04 


-0 


.018 


± 





.002 





.50 


± 


0. 


.06 


0.; 


3 ± 0.1 


1.3 ± 0.1 


800. 


-0 


.38 


± 


0. 


.04 


-0 


.023 


± 


0. 


.002 





.60 


± 


0. 


.07 


0.9 ± 0.1 


1.5 ± 0.1 


1000. 


-0 


.40 


± 


0. 


.02 


-0 


.024 


± 


0. 


.001 





.61 


± 


0. 


.03 


0.93 


± 0.05 


1.53 


± 0.08 


1100. 


-0 


.43 


± 


0. 


.03 


-0 


.025 


± 


0. 


.001 





.63 


± 


0. 


.04 


0.99 


± 0.06 


1.57 


± 0.10 


1500. 


-0 


.45 


± 


0. 


.02 


-0 


.026 


± 


0. 


.001 





,64 


± 


0. 


.02 


0.97 


± 0.04 


1.58 


± 0.06 


1700. 


-0 


.46 


± 


0. 


.02 


-0 


.027 


± 


0. 


.001 





.66 


± 


0. 


.03 


0.99 


± 0.05 


1.61 


± 0.08 


2000. 


-0 


.49 


± 


0. 


.01 


-0 


.028 


± 


0. 


.001 





.68 


± 


0. 


.02 


1.03 


± 0.03 


1.66 


± 0.05 


3000. 


-0 


.55 


± 


0. 


.02 


-0 


.031 


± 


0. 


.001 





.72 


± 


0. 


.02 


1.07 


± 0.03 


1.71 


± 0.05 


4000. 


-0 


.60 


± 


0. 


.02 


-0 


.033 


± 


0. 


.001 





.76 


± 


0. 


.02 


1.12 


± 0.03 


1.77 


± 0.04 


5000. 


-0 


.61 


± 


0. 


.02 


-0 


.033 


± 


0. 


.001 





.76 


± 


0. 


.02 


1.11 


± 0.03 


1.75 


± 0.04 


6000. 


-0 


.65 


± 


0. 


.01 


-0 


.035 


± 


0. 


.001 





,89 


± 


0. 


.01 


1.17 


± 0.02 


1.82 


± 0.02 


7000. 


-0 


.67 


± 


0. 


.02 


-0 


.035 


± 


0. 


.001 





.81 


± 


0. 


.02 


1.18 


± 0.03 


1.84 


± 0.04 


8000. 


-0 


.69 


± 


0. 


.01 


-0 


.036 


± 


0. 


.001 





.83 


± 


0. 


.01 


1.21 


± 0.02 


1.88 


± 0.03 


9000. 


-0 


.70 


± 


0. 


.02 


-0 


.036 


± 


0. 


.001 





.83 


± 


0. 


.01 


1.22 


± 0.02 


1.90 


± 0.03 


10000. 


-0 


.70 


± 


0. 


.01 


-0 


.036 


± 


0. 


.001 





.84 


± 


0. 


.01 


1.23 


± 0.02 


1.92 


± 0.03 


11000. 


-0 


.69 


± 


0. 


.01 


-0 


.035 


± 


0. 


.001 





,84 


± 


0. 


.02 


1.23 


± 0.01 


1.91 


± 0.02 
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-0 


.69 


± 


0. 


.01 


-0 


.035 


± 


0. 


.001 





.85 


± 


0. 


.01 
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± 0.02 


1.94 


± 0.03 
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-0 


.71 


± 


0. 


.02 


-0 


.036 


± 


0. 


.001 





.86 


± 


0. 


.01 


1.27 


± 0.02 


1.97 


± 0.03 


14000. 


-0 


.72 


± 


0. 


.01 


-0 


.036 


± 


0. 


.001 





.87 


± 


0. 


.01 


1.28 


± 0.02 


1.98 


± 0.03 


15000. 


-0 


.72 


± 


0. 


.02 


-0 


.036 


± 


0. 


.001 





.87 


± 


0. 


.02 


1.28 


± 0.02 


1.97 


± 0.04 



with age, also because a larger number of stars is needed 
to reach the given total magnitude. 

4.2. IMF variations 

To evaluate the effect of IMF variations on integrated 
colours we changed the Kroupa (2002) IMF exponent (x) 
within the estimated uncertainty for masses M > 0.5 Mq. 
In Fig. O models with My 1 — —8 calculated with the 
lower, central and upper value of the IMF exponent are 
plotted in selected pass-bands; the corresponding data ta- 
bles are available as electronic tables at CDS. 

As obtained for metal-rich models investigated in pre- 
vious papers Csee e.g.lMarastonlll998UBrocato et"aTll2000l 
lYil l2003t iBruzual fc Chariot! l2003|k we find that colours 
variations are negligible also for very low-metallicity mod- 
els, at least in the range of exponent investigated here, 
being the colours variations well within the intrinsic un- 
certainty, for both values of the brightness of the cluster 
(My 1 = —4 mag and —8). The errorbars behavior is sim- 
ilar to what described above: colour fluctuations increase 
from optical-optical to optical-NIR colours and for ages 
^ 1 Gyr. 

Changing the distribution of very low mass stars 
(M < 0.5 Mq) has negligible effects on photometric in- 
dices whatever the total visual magnitude, as already 



shown by iBrocato et al.l <)2000|) for ages greater than 5 
Gyr. In poorly populated clusters (My* — —4 mag) as in 
the richest ones (My 1 = — 8 mag), reducing the number 
of stars with M < 0.5 Mq does not affect the integrated 
colours since the contribution of these sta rs to the total 
V lig ht is only of the order of few percent ijBrocato et alJ 
2000). The contribution of such stars slightly increases 
with age, but it remains comparable to the statistical er- 
rors even at the oldest ages (at least in the range explored 
here) . 

4.3. Comparison with previous works 

Table 0] summarizes the main ingredients used by various 
authors in their synthesis code, and Fig. compares the 
time-evolution of B — V and V — K colours of the present 
work with similar recent results available in the literature. 
In the figure our models are plotted with la intrinsic error, 
while other autho rs do not estimate th e colours statistical 
errors, except for lBrocato et alJ l|2000f) . If colours are not 
available for exactly the same metallicity value adopted 
in the present work (Z = 0.0002), the two closest Z val- 
ues are plotted. The models referred as Anders & Fritze-v. 
Alvensleben (2003) are the ones of the quoted paper, up- 
dated co ncerning emission l ines by private communication 
(see also lSchulz et al.ll2002|) . 
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Fig. 3. Time evolution of selected integrated colours in the standard UBVRIJHK filters for models with My* = —4 
(red triangles, dashed lines) and My* = —8 (black points, solid lines). For each model 1 a colour dispersion is reported, 
see text. 



The present pre dictions well agree with those by 
iBrocato et all ll2000l). There is also a general good agree- 
ment with iKurth et all l|l999|) . and iBruzual fe Chariol 
l)2003(l . both based on Padua 1994 stellar evolutionary 
tracks (see the quoted papers for the detailed list of 
references). Differences can be found in the age range 
8.6 ^ log age(yr) ^ 9 for optical-NIR colours, possibly 
due to the tr e atmen t of the AGB and TP-AGB phases. 
iGirardi et al.l l)2000|) and Anders & Fritze-v. Alvensleben 
(2003) predict redder colo urs than all the othe r authors for 
age > 100 Myr, as well as lZhang et"aD ll2002l) for age Z 3 
Gyr. Differently from the others, both iMarastonl (|2005|) 
and our models are based on no-overshooting tracks. 
Nevertheless, while the present B — V colours well agree 
with Maraston's models at all ages in common, V — K 
predictions are bluer than those by IMarastonl l)2005[) at 
ages younger than ~ 3 Gyr. 

Fig.|S|shows that the B — V and V—K predictions from 
different authors agree within 0.1 mag at log age ^ 9.3. 
Large model-to-model differences arise at young and inter- 
mediate ages, when the contribution of AGB and TP-AGB 
stars to the total flux is high, especially in the NIR bands 
(see for a discussion IMarastonl Il99st IGirardi k, Bertellil 
1998). This has a direct consequence on predicted colours, 



since a large variety of stellar ingredients and prescriptions 
are used to simulate these phases. Temperature-colours re- 
lations, mass-loss prescriptions, and analytical description 
of the TP evolution all affect the photometric properties 
of such kind of stars, thus it is not straightforward to 
individuate a single origin (or multiple ones) to explain 
differences shown in Fig. in the intermediate-age range. 

For instance, mass-loss processes affecting TP-AGB 
stars are one of the physical mechanisms largely affect- 
ing the TP-AGB stars evolution. However, the effect of 
mass loss on the observational properties of TP stars are 
largely unknown to date. To give an indication on how 
mass-loss may affects colour predictions, we made numer- 
ical experiments by computed integrated colours in the 
extreme case when no TP-AGB stars are present in the 
population, to mi mic in such a way a very efficient mass- 
loss rate (see also IMarastonl Il998). As expected, the ef- 
fect is larger for intermediate- age populations, and in the 
optical-NIR colours. The V — K colour becomes bluer than 
the standard value by ~ 0.5 mag at age = 500 Myr, 
and the B — V colour decreases only by 0.06 mag at the 
same age. This effect tends to vanish by increasing the 
age. On the other hand, if we lower the mass-loss effi- 
ciency with respect to our standard assumption (BH) by 
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Fig. 5. Time evolution of select e d inte grated colours for SSPs with My' = —8 calculated with the lower, central 
and upper value for the iKrouoal l|2002|) IMF exponent (a) for M > 0.5M Q . Black dots indicate our reference model 
(a = 2.3), red triangles indicate models with a — 2.0 and blue stars models with a = 2.6. 



using a Reimers' law with r\ of the order of 1, and V — K 
colours redder than our 'standard' models are obtained. 
As a consequence the colour decline at age corresponding 
to the appearance of He-core degenerate stars is more evi- 
dent i t ~ 1 Gy r\ similarly to what ob tained bv lMarastonl 
l)2005|) . and bv lRaimondo et all l)2005|) for surface bright- 
ness fluctuation colours. 



5. Mass-to-light ratio 

Table [H] lists our predictions for the mass-to-light ratio 
{M. I L) in various photometric bands as function of age for 
Mtf 1 = -8; similar tables for M#" = -6 and M^ ot = -4 
are available at CDS. 

The total mass M. includes all the evolving stars up 
to the AGB tip (M > 0.1 M Q ) and WDs. Stars with mass 



higher than M up = 7M Q are assumed to leave a neu- 
tron star (NS) as a remnant. To evaluate the maximum 
mass fraction of NS we set th e neutron star mass at its 
upper limit (« 2M Q , see e.g. iBombaci et aDl20D4h . and 
we choose high cluster age (« 12 Gyr); we found that 
the NS mass fraction constitutes at maximum few per- 
cent of the cluster mass in agreem ent with the calcula- 
tions bv IVesoerini fc Heggid <)l997l) who found an upper 
value of w 1%. The contribution of massive black holes to 
the total mass depe nds on assumpti on on the mass of the 
remnant and IMF l)Marastonll2005|). By using the s ame 
prescriptions on black holes mass as lMarasto we 



estimated that for a Kroupa IMF the mass fraction of both 
neutron stars and black holes does not exceed 5-10%. 

In this paper, the Ai/L values are calculated with- 
out taking into account any dynamical processes, e.g. 
evaporation of st ars due to two body relaxation or disk 
shocking (see e.g . ISpitzerlfl987t IVesperini fc Heggid lT997t 
iBoilv et al.ll2005fl . A treatment of these phenomena is be- 
yond the scope of this work; N-body simulations show- 
ing the effects of this processes on the total mass and 
IMF shape as a function of time, and on the galac- 
tocentric distance of the cluster can be found in e.g. 
IVesperini fc Hegglel (|l99lt . We only note that since stel- 
lar evaporation is more efficient for low mass stars, one 
expects a flattening of the IMF as the dynamical evolu- 
tion of a cluster proceeds. Evaluations of the mass fraction 
of WDs expected to be lost fr om clusters due to dynami- 
cal evap oration can be found i nlYesi^firinL^JIeggie] ( 199^ 
(see alsolFellhrnier et aJ2003tlHurlev fc Sba,ral200.^. The 
presence of binary stars is not taken into account as well; 
the lack of information about the binary frequency, the 
distribution of binaries of different mass ratios, the sepa- 
ration of the components, and the occurrence of explosion 
of novae and supernovae prevents a quantitative treat- 
ment of this phenomenon. An attempt to include binary 
systems in p opulation syn t hesis models has recently been 
presented bv lZhang et alJ l)2005j) . 

In Fig. we show the influence of VLM stars (M < 
0.6 Mq) and WDs on the Ai/L, L is the bolometric lumi- 
nosity in solar units, at fixed absolute visual magnitude 
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Table 4. The physical ingredients adopted in the models plotted in Fig|S| The quoted values for the overshooting 
(OS) efficiency, expressed in unit of pressure-scale-height above the Schwarzschild convective border, are only indicative 
because different sets of evolutionary tracks adopt different prescriptions for the overshooting mechanisms and for the 
minimum mass at which the overshooting is fully efficient (see the related papers for more details); Y=YES; N=NO. 
WD indicates if white dwarfs are included in the calculations; Ai/L if mass-to-light ratio values are available. In 
the paper by Kurth et al. (1999) the normalization procedure is not described. In the paper by Anders & Fritze-v. 
Alvensleben (2003), the SSPs have been evolved to the selected ages starting with an initial stellar mass of 10 6 Mq; 
the adopted value for M up is not specified. 1 Detailed information on the adopted evolutionary tracks can be found 
in the paper quoted in Col. 1. 2 Latest version available on the web, see Maraston (2005). 3 Sinthetic models are also 
calculated for Salpeter (1955) IMF. 4 See discussion in Maraston (2005). 5 Carbon stars included. 
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(My = —8). As already known, VLM and WDs provide a 
nearly negli gible contributi on to the total luminosity (see 
for example lMarastonll998|) . while give a significative con- 
tribution to the total mass. Fig.0shows that the fraction 
of low mass star s and white dwa rfs increases with the clus- 
ter age (see also lMarastonll998| . Note that, since the sim- 
ulations are performed at fixed absolute visual magnitude, 
the total mass of the cluster varies with age, as shown in 
Fig. 03 where we plot the case of My — —8. This is be- 
cause as the age increases the mass of the typical star in 
the cluster decreases. 



Changing M up affects the number of stars cooling as 
WDs, or in other words the mass fraction locked as WDs, 
and thus the cluster total mass. Different input ingredi- 
ents of stellar evolutionary tracks, as e.g. the neutrino 
energy losses in the core, lead to different predictions 
of My. n at fixed metalli city, see e.g. iTornambe fc Chieffil 
HmL M..„ = 5.7M^VlPols et al.l l|l998i M, m = 5Mm). 
iDominguez et all ()l999t M up = 6.5 Mq), ICariulo et alJ 
l|2004 M„„ = 5 4- 6 Mq). We explore the effect of chang- 
ing M up . We find that increasing M up from 5 up to 7 Mq 
does not cause variations in the mass-to-light ratios, at 
least for the adopted IMF. Even pushing M up up to IOMq 
does not cause sensitive change in the mass-to-light ratios 
(Table 0). 
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Fig. 8. Time behaviour of the logarithm of the cluster 
total mass needed to reach My * = —8 normalized to the 
maximum value, as a function of age. 

5.1. My 1 sensitivity 

Fig. El shows the time evolution of M/L ratio in selected 
passbands for standard models with My 1 = —8 (Table 
0) and —4 mag. The corresponding data tables for other 
passbands and for the case of My 1 — — 6 are available at 
CDS. Note that the mass-to-light ratios generally increase 
with age independently of the photometric band, because 
the total mass increases. This result is similar to what ob- 
tained by other authors, who do not keep constant the lu- 



Fagiolini et al.: Monte Carlo Simulations of Metal Poor Star Clusters 



13 



■ Kurth ei al. 1999, Z=0.0004 
. Kurth et al. 1999, Z=0.0001 
- Bruzual & Chariot 2003, Z-0.0001 

• Bruzual & Chariot 2003, Z-0.0004 

• Present work, Z- 0.0 002 




3.5 9 9.5 

log[age(yr)] 




Kurth et al. 1999, Z=0.0001 
Kurth et al. 1999, Z=O.0OD4 
Bruzual & Chariot 2003, Z=0.0001 
Bruzual & Chariot 2003, Z=0.0004 
Present work, Z=0.0002 



1. Mlaraston (2005) 

• Zhang et al. 2002, Z-0.0002 

Anders S Fritze-v Alvensleben 2003, Z=0.0004 
- -■ Girardi et al. 2000, Z=0.0004 
Brocato et al. (2000) Z=0.0001 

• Present work, Z=0.0002 




8.5 9 9.5 

log[age(yr)] 




a Maraslon (2005) 

• Zhang et al. 2002, Z=0.0002 

Anders & Fritze-v Alvensleben 2003. Z=0.0004 
. Girardi et al. 2000, Z=0.0004 

• Present werk, Z=0. 0002 



9.5 



10 



Fig. 6. Comparison between theo retical integrated co lours (B — V, upper panels; V — K, lower panels) f rom recent pa- 
pers: p resent mod els (black dots1:lKurth et alJ lIlQQfll red Z=0. 0001. green Z =0.0004 squares, solid line). lBrocato et all 
cyan dots). lAnders fc Fritze-v. Alvensleben I <)2003l orange solid line), C-irardi et al.1 l)2000l blue squares, short- 



dashed linel. lZhang et alj (I2Q02L ' pink circles, dotted line'l. lBruzual fc Charloll pOOa violet Z=0.0001. brown Z=0.0004 
circles, dashed linel. lMarastonl (f2005l red triangles, dashed line). 
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Fig. 7. Left panel: Fraction of very low mass stars (M < 0.6 Mq) that contribute to the mass-luminosity ratio (A4/L) 
(in solar units) as a function of the age for My * = —8. Filled circles represent our standard model in which all masses 
with M > 0.08 Mq are included; the effect of considering only stars with M > 0.4M Q (filled triangles) or M > 0.6 Mq 
(filled squares) is analyzed separately. Right panel: as in the left panel, but when the WD population is (filled circles, 
standard model) or is not (filled squares) included. 



minosity Ly (see e.g. iBruzual fc Charlotll2003l iMarastoi] 
HH). 

The only effect of changing the absolute magnitude 
from My = —8 to —4 on mass-to-light ratios in the optical 
bands is to increase the intrinsic uncertainties. On the 
contrary, other than this effect mass-to-light ratios in the 
NIR bands suffer a larger scatter due to the decrease of red 
giants contributors to the total NIR luminosity, as shown 
in the right lower panel of the figure for if -band. 



Since is fixed, the behaviour of (Ai/Ly) as a 

function of the age reflects the trend of the cluster mass 
(Fig. |SJ); it is quite linear for ages > 1 Gyr in all pass- 
bands indicating it is driven by the grown of the total 
mass. For ages < 1 Gyr NIR colours are flatter than the 
optical (U,B,V) ones due to the increase of the infrared 
luminosity following the develop of an extended asymp- 
totic giant branch. We also find a dimming in M/Lk at 
the age corresponding to the AGB phase-transition con- 
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Table 5. Mass- to -light ratios for selected ages in several photometric passbands from our standard model with different 
assumptions of M up . 



M/Ly M/Lu M/ Lb M/Lr M/Li M/Lj M/L K 

M U p = 5 M Q 

500 Myr 0.21 ± 0.01 0.126 ±0.003 0.139 ±0.004 0.25 ± 0.02 0.28 ± 0.04 0.28 ± 0.06 0.26 ± 0.08 

800 Myr 0.28 ± 0.02 0.194 ±0.007 0.210 ± 0.008 0.31 ± 0.03 0.31 ± 0.04 0.29 ± 0.05 0.25 ± 0.06 

3Gyr 0.66 ±0.03 0.55 ± 0.01 0.58 ± 0.02 0.67 ± 0.03 0.64 ± 0.03 0.58 ± 0.04 0.50 ± 0.04 

WGyr 1.55±0.04 1.45 ±0.02 1.58 ±0.03 1.47 ±0.04 1.33 ±0.05 1.15 ±0.05 0.95 ± 0.05 
M up = 6.5 M Q 

500 Myr 0.22 ± 0.01 0.127 ± 0.003 0.140 ± 0.004 0.26 ± 0.02 0.28 ± 0.03 0.29 ± 0.06 0.26 ± 0.08 

800 Myr 0.28 ± 0.01 0.191 ± 0.005 0.207 ± 0.006 0.30 ± 0.02 0.30 ± 0.03 0.28 ± 0.05 0.25 ± 0.06 

3Gyr 0.66 ± 0.03 0.55 ± 0.01 0.58 ± 0.02 0.67 ± 0.03 0.64 ± 0.04 0.58 ± 0.04 0.49 ± 0.04 

WGyr 1.58 ± 0.05 1.47 ± 0.02 1.60 ± 0.03 1.49 ± 0.05 1.36 ± 0.06 1.17 ± 0.06 0.97 ± 0.06 

M up = 10 M Q 

500 Myr 0.22 ± 0.01 0.130 ±0.002 0.143 ±0.003 0.26 ± 0.02 0.28 ± 0.04 0.28 ± 0.06 0.26 ± 0.08 

800 Myr 0.29 ± 0.01 0.199 ±0.005 0.214 ± 0.005 0.31 ± 0.02 0.31 ± 0.03 0.29 ± 0.05 0.25 ± 0.05 

3Gyr 0.68 ±0.03 0.56 ± 0.01 0.60 ± 0.02 0.69 ± 0.04 0.66 ± 0.04 0.60 ± 0.04 0.51 ± 0.05 

WGyr 1.58 ±0.05 1.48 ±0.02 1.61 ±0.04 1.49 ± 0.06 1.35 ±0.06 1.16 ±0.06 0.96 ± 0.06 
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Fig. 9. Time behaviour of the Ai/L ratio (in solar units) for selected passbands (U, B,V, K). For each model la 
dispersion is shown. Red triangles (dashed line) indicate models with My 1 — —4, and black dots (continuos line) 
models with M\f - 8. 



firming the result found by iMarastonl J2005) for higher very satisfactory, while lAnders k. Fritze-v.Alvensleben I 

, 11- •, Innnri, l- , T l ,1 1 • 1 1 , 1 o r- 



metallicity 

Finally, 
ues with 



Fig. ma 

the ones 



co mpares the 
bv iBruzual 



prese nt val- 



lAnders fc Fritze-v.Alvensleben I <|200^ . and iMarastorJ 
(2005) in two photometric bands, namely B and 
K. Concerning t h e .B- band, the agree m ent with 
IBruzual fc Charlod (|2003|) and IMarastonl l|2005|) is 



1200.^ predict slightly higher values at log age ^ 8.5. 



A more complex behaviour is shown by JVI/Lk- at 
log age <; 9.3 the A4/Lk is independent of the model, 
while model-dependent at younger ages. This reflects the 
great uncertainties in simulating the AGB phase we have 
already discussed for colours. 
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Fig. 10. Time evolution of Ml Lb ( left panel) and A4/Lk (right panel) in our reference model (black dots) compared 
with the iBruzual fc Chariot! l|2003l B&C 2003) for models with Z = 0.0001 (violet dashed line), and Z = 0.0004 
(brow n dashed line), Anders & Fritze-v. Alvensleben (2003, A&F 2003) models (orange solid line), and iMarastonl 
(2005, M2005) for Z = 0.0003 (red triangles, dashed line). The IMF adopted by various authors are reported in 
Table H 

Table 6. Theoretical Ai/L as a function of the age in different photometric bands (solar units) obtained for standard 
models with M\f = -8. 
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5.2. IMF variations 

To evaluate the effect of the IMF shape on Ai/L we 
changed the IMF exponents within the estimated uncer- 



tainty l|Kroupal2002() . Fig.[TT]and Fig.dshow the effect of 
changing the IMF exponents with respect to the Kroupa's 
formulation for masses higher and lower than 0.5 M Q , re- 
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Fig. 11. Time evolution of M./L (in solar units) in sel ected passband s (U, B,V, K) for models with My 1 = —8 
calculated with the lower, central and upper value for the iKrouoal l|2002j) IMF exponent for M < 0.5M©. Black dots 
indicate our reference model (a ~ 1.3), red triangles indicate models with x=0.8 and blue stars models with a = 1.8. 



spectively. Models with My* = —8 mag are plotted in 
selected passbands while results for other passbands and 
for My' = — 4 mag are available at CDS. 

The A4/L ratio is more sensitive than integrated 
colour s to the IMF shape, as alr e ady noted bv iMarastonl 
l|l998|) : lBruzual fc Charloll|2003|) : lMarastonll)2005|) . From 
Fig. El it is evident that, by decreasing the IMF slope at 
M < 0.5M@, the mass of the system decreases, while 
the luminosity remains almost constant. The final effect is 
to have smaller M/L ratios than in the standard models 
(a = 1.3). At age lower than few billion years, the un- 
certainty due to stochastic effects is comparable or larger 
than the effects of exponent variations. Indeed, the IMF 
variations can be appreciated at ages of the order of 
10 10 yr. We also find that the uncertainty predicted for 
My' = —4 models prevents to detect significant varia- 
tions at any age. 

A different situation is found in Fig. ^] The mass-to- 
light ratios slightly depend on IMF for age younger than 
a few billion years, while they are nearly insensitive at 
old age. This is because an IMF flatter than a — 2.3 at 
M > 0.5 (red triangles) requires a total mass smaller than 
in the standard case (a = 2.3, black circles) to reach the 
given luminosity, thus predicting low M./Ly. The effect is 
enhanced in young stellar populations in which the total 
mass is dominated by massive stars. On the other hand, by 
increasing the exponent up to 2.6 (blue stars), the cluster 
mass needed to have My = — 8 increases, leading to high 



M./L values. Again, mass-to-light ratios in the NIR bands 
are more affected by stochastic phenomena. 

Finally, note that in the case of small variations (a = 
2 -T- 2.6) explored here, the contribution of massive rem- 
nants to the total mass does not significantly change, while 
it is shown to be effective if large variations (a = .5-^3.5) 
are adopted, see for a discussion M arastonl l|1998|) . 

6. Summary and conclusions 

In this work we have analyzed the intrinsic uncertainties 
due to stochastic effects on integrated colours and mass- 
to-light ratios of metal-poor stellar clusters as a function 
of the total visual magnitude. The calculations are per- 
formed for three different values of My * and for a fine 
grid of stellar ages. Statistical errors are shown to be cru- 
cial, especially in the extreme case My' = —4 they are so 
high as to prevent precise quantitative evaluations for all 
ages. Calculations are done in the standard UBVRIJHK 
photometric passbands and in the Hubble Space Telescope 
bands (WFPC2 and NICMOS systems). 

We checked the consistency of our models on observa- 
tional properties of three metal poor clusters, namely M68, 
M15, and M30, which differ mainly in the absolute visual 
magnitude, and HB morphology. For each cluster we were 
able to reproduce both the features of the observed CMD 
and the integrated colours, showing that the HB morphol- 
ogy does not influence the photometric indices taken into 
account. 
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Fig. 12. Time evolution of M./L (in solar units) in sel ected passband s (U, B,V, K) for models with My 1 = —8 
calculated with the lower, central and upper value for the lKroupal l|2002() IMF exponent for M > 0.5M Q . Black dots 
indicate our reference model (a = 2.3), red triangles indicate models with a = 2.0, and blue stars indicate models 
with a — 2.6. 



A comparisons with recent results available in the liter- 
ature shows, in some cases, non-negligible differences due 
to the large variety of prescriptions used in the model cal- 
culations. 

The uncertainties in the results on both colours and 
mass-to-light ratios due to the still present uncertainty on 
the IMF shape have been quantitatively estimated; while 
the colours fluctuations remains within theoretical uncer- 
tainties the Ai /L ratio is more sensitive to the IMF shape. 
We also showed that the influence on M. / L of the adopted 
value for the minimum mass for which carbon burning is 
ignite is quite negligible. 

Acknowledgements. This work is dedicated to the memory of 
Vittorio Castellani, whose advise always supported us, and 
whose enthusiasm and dedication in challenging the astro- 
physics questions are an invaluable example for all of us. 
We warmly thank E. Brocato, P.G. Prada Moroni and S.N. 
Shore for useful discussions and for a careful reading of the 
manuscript. We are grateful to the anonymous referee for 
her/his suggestions and comments that greatly improved the 
paper. Financial support for this work was provided by MIUR- 
COFIN 2003. This work made use of computational resources 
granted by the Consorzio di Ricerca del Gran Sasso according 
to the Progetto 6 'Calcolo Evoluto e sue Applicazioni (RSV6)'- 
Cluster Cll/B. This paper utilizes the HST-snapshot database 
by the Globular Cluster Group of the Padova Astronomy 
Department, and the Catalog of parameters for Milky Way 
Globular Clusters by W.E. Harris. 



References 

Allard F., Hauschildt P.H., Alexander D.R., Starrfield S., 

1997, ARA&A, 35, 137 
Alves D.R., Bond H.E., Livio M., 2000, AJ, 120, 2044 
Anders P., Fritze-v. Alvensleben U., 2003, A&A, 401, 1063 
Angeletti L., Dolcetta R., Giannone P., 1980, Ap&SS, 69, 

45 

Bailyn, CD., Grindlay, J.E., Cohn, H., Lugger, P.M., 
Stetson, P.B., Hesser J.E. 1989, AJ, 98, 882 

Baraffe I., Chabrier G., Allard F., Hauschildt P.H., 1997, 
A&A, 327, 1054 

Barbara C, Bertelli C, 1977, A&A, 54, 243 

Barmina R., Girardi L., Chiosi C, 2002, A&A, 385, 847 

Baud B., Habing H.J., 1983, A&A, 127, 73 

Bergeron P., Saumon D., Wesemael F., 1995, ApJ, 443, 
764 

Bessell M.S., Castelli F., Plez B., 1998, A&A, 333, 231 
Blocker T., Herwig F., Driebe T., Bramkamp H., 

Schonberner D., 1997, In: IAU Symp. 180: Planetary 

Nebulae, 389 

Boily CM., Lancon A., Deiters S., Heggie D.C, 2005, 

ApJ, 620, L27 
Bombaci I., Parenti I., Vidaha I., 2004, ApJ, 614, 314 
Bressan A., Chiosi C, Fagotto F., 1994, ApJS, 94, 63 
Brocato E., Caputo F., di Giorgio A.M., Santolamazza P., 

Richichi A., 1990, Memorie della Societa Astronomica 

Italiana, 61, 137 



18 



Fagiolini ot al.: Monte Carlo Simulations of Metal Poor Star Clusters 



Brocato E., Matteucci F., Mazzitelli I., Tornambe A., 

1990b, ApJ, 349, 458 
Brocato E., Castcllani V., Raimondo G., Romaniello M., 

1999, A&AS, 136, 65 

Brocato E., Castellani V., Poli F.M., Raimondo G., 2000, 

A&AS, 146, 91 
Brocato E., Castcllani V., Di Carlo E., Raimondo G., 

Walker A.R., 2003, AJ, 125, 3111 
Bruzual G., Chariot S., 2003, MNRAS, 344, 1000 
Burstein D., Faber S. M., Gaskell C. M., Krumm N., 1984, 

ApJ, 287,586 
Buzzoni A., 1993, A&A, 275, 433 

Cardclli J.A., Clayton G.C., Mathis J.S., 1989, ApJ, 345, 
245 

Cariulo P., Degl'Innocenti S., Castellani V., 2004, A&A, 
421, 1121 

Carretta E., Gratton R.G., Clementini G., Fusi Pecci F., 

2000, ApJ, 533, 215 

Castellani V., Chieffi A., Tornambe A., Pulone L., 1985, 
ApJ, 296, 204 

Castellani V., Degl'Innocenti S., Marconi M., Prada 

Moroni P.G., Scstito P., 2003, A&A404, 645 
Castelli F., 1999, A&A, 346, 564 

Castelli F., Gratton R.G., Kurucz R.L., 1997, A&A, 318, 
841 

Cervino, M. and Valls-Gabaud, D., 2003, MNRAS, 338, 
481 

Chabrier G., Baraffe I., 1997, A&A, 327, 1039 
Chabrier G., Mera D., 1997, A&A, 328, 83 
Chariot S., Bruzual G., 1991, ApJ, 367, 126 
Chariot S., Worthey G., Bressan A., 1996, ApJ, 457, 625 
Chiosi C, Bertelli G., Bressan A., 1988, A&A, 196, 84 
Ciacio F., Degl'Innocenti S., Ricci B., 1997, A&AS 123, 
449 

Cote P., Marzkc R.O., West M.J., 1998, ApJ, 501, 554 
De Marchi G., Paresce F. 1994, ApJ, 422, 597 
Di Criscienzo M., Marconi M., Caputo F., 2004, ApJ, 612, 
1092 

Dominguez I., Chieffi A., Limongi M., Straniero O., 1999, 
ApJ, 524, 226 

Fagiolini M., 2004, Diploma Thesis, University of Pisa, 
Italy 

Fellhauer M., Lin D.N.C., Bolte M., Aarseth S.J., 

Williams K.A., 2003, ApJ, 595, L53 
Ferraro F.R., Messineo M., Fusi Pecci F., et al., 1999, AJ, 

118, 1738 

Girardi L., Bertelli L., 1998, MNRAS,300, 533 
Girardi L., Bressan A., Bertelli G., Chiosi C, 2000, A&AS, 
141, 371 

Harris W.E., 1996, VizieR Online Data Catalog, 7195, 
Harris W.E., 2003, In: A Decade of Hubble Space 

Telescope Science, 78-100 
Hurley J.R., Shara M.M., 2003, ApJ, 589, 179 
Kron G.E., Mayall N.U., 1960, AJ, 65, 581 
Kroupa P., 2002, Science, 295, 82 

Kurth O.M., Fritzc-v. Alvcnslebcn IT., Fricke K.J., 1999, 

A&AS, 138, 19 
Lattanzio J. C, Wood P. R., 2003, in Asymptotic giant 



branch stars, by Harm J. Habing and Hans Olofsson. 

Astronomy and astrophysics library, New York, Berlin: 

Springer, 2003, p. 23 (2003) 
Larsen S.S., 2000, MNRAS, 319, 893 
Lee Y., Demarque P., Zinn R., 1994, ApJ, 423, 248 
Maeder A., Zahn J., 1998, A&A, 334, 1000 
Maraston C, 1998, MNRAS, 300, 872 
Maraston C, 2005, MNRAS, 362, 799 
Marigo, P., Girardi, L., & Bressan, A. 1999, A&A, 344, 

123 

Matteucci A., Ripepi V., Brocato E., Castellani V., 2002, 

A&A, 387, 861 
McNamara B.J., Harrison T.E., Baumgardt H., 2004, 

ApJ, 602, 264 
Origlia L., Leitherer C, 2000, AJ, 119, 2018 
Palacios A., Talon S., Charbonnel C, Forestini M., 2003, 

A&A, 399, 603 
Pasquali A., De Marchi G., Pulone L., Brigas M. S., 2004, 

A&A, 428, 469 
Pease F.G., 1928, PASP, 40, 342 

Piotto G., King I.R., Djorgovski S.G., et al., 2002, A&A, 
391, 945 

Pols O.R., Schroder K., Hurley J.R., Tout C.A., Egglcton 
P.P., 1998, MNRAS, 298, 525 

Prada Moroni P.G., Straniero O., 2005, ApJ, in publica- 
tion 

Pryor C, Meylan G., 1993, In: ASP Conf. Ser. 50: 
Structure and Dynamics of Globular Clusters, 357 

Raimondo G., Brocato E., Cantiello M., Capaccioli M., 
2005, AJ, 130, 2625 

Reed B.C., 1985, PASP, 97, 120 

Reed B.C., Hcsser J.E., Shawl S.J., 1988, PASP, 100, 545 
Rcimcrs D., 1975, Memoires of the Societe Royale des 

Sciences de Liege, 8, 369 
Rcnzini, A., & Buzzoni, A. 1986, in Proceedings of the 

Fourth Workshop, Erice, Italy, March 12-22, 1985, 

Spectral evolution of galaxies (Dordrecht, D. Rcidcl 

Publishing Co.), 195 
Renzini A., Voli M., 1981, A&A, 94, 175 
Rey S., Yoon S., Lee Y., Chaboyer B., Sarajedini A., 2001, 

AJ, 122, 3219 

Richard O., Michaud G., Richer J., 2002, ApJ, 580, 1100 
Richer J., Michaud G., Rogers F., et al., 1998, ApJ, 492, 
833 

Rosenberg A., Aparicio A., Saviane I., et al., 2000, A&AS, 
145, 451 

Salaris M., Garcfa-Berro E., Hernanz M., Isern J., Saumon 

D., 2000, ApJ, 544, 1036 
Sandquist E.L., Bolte M., Langer G.E., Hesser J.E., 

Mendes de Oliveira C, 1999, ApJ, 518, 262 
Santos J.F.C., Frogel J. A., 1997, ApJ, 479, 764 
Saumon D., Jacobson S.B., 1999, ApJ, 511, L107 
Schulz J., Fritze-v. Alvenslcben U., Moller C.S., Fricke 

K.J., 2002, A&A, 392, 1 
Spitzcr L., 1987, Dynamical evolution of globular clusters, 

Princeton, NJ, Princeton University Press, 1987, 191 p. 
Stetson P.B., 1994, PASP, 106, 250 

Straniero O., Domngucz, I., Cristallo, R., Gallino, R. 2003, 



Fagiolini et al.: Monte Carlo Simulations of Metal Poor Star Clusters 



PASA, 20, 389 
Thoul A.A., Bahcall J.N., Loeb A., 1994, ApJ, 421, 828 
Tinsley B.M., 1972, A&A, 20, 383 
Tornambe A., Chieffi A., 1986, MNRAS, 220, 529 
Vazdekis A., 1999, ApJ, 513, 224 
Vesperini E., Heggie D.C., 1997, MNRAS, 289, 898 
Wagenhuber J., Groenewegen M.A.T., 1998, A&A, 340, 

183 

Walker A.R., 1994, AJ, 108, 555 

West M.J., Cote P., Marzke R.O., Jordan A., 2004, 

Nature, 427, 31 
Worthcy G., 1994, ApJS, 95, 106 
Yi S.K., 2003, ApJ, 582, 202 

Zhang F., Han Z., Li L., Hurley J.R., 2002, MNRAS, 334, 
883 

Zhang F., Li L., Han Z. 2005, MNRAS, 364, 503 



